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Abstract 

We establish that a mode-coupling approximation for the dynamics of multi-component systems 
obeying Smoluchowski dynamics preserves a subtle yet fundamental property: the matrices of 
partial density correlation functions are completely monotone, i.e. they can exactly be written 
as superpositions of decaying exponentials only. This statement holds, no matter what further 
approximations are needed to calculate the theory's coupling parameters. The long-time limit of 
these functions fulfills a maximum property, and an iteration scheme for its numerical determination 
is given. We also show the existence of a unique solution to the equations of motion for which 
power series both for short times and small frequencies exist, the latter except at special points 
where ergodic-to-nonergodic transitions occur. These transitions are bifurcations that are proven 
to be of the cuspoid family. 

PACS numbers: 82.70.Dd, 05.40.-a, 02.10.Sp 
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I. INTRODUCTION 



Density correlation functions are convenient tools to characterize the dynamics of liquids 
or disordered systems. They can be measured in experiment, e.g. by inelastic neutron 
scattering, dynamic light scattering in colloidal systems, or can be determined from computer 
simluation techniques. On the other hand, one can calculate them from theory, but in the 
case of strongly interacting systems, one typically has to invoke certain approximations in 
order to obtain the desired results. There are, however, some general properties of such 
correlation functions, directly related to the time-evolution operator of the system. It is a 
nontrivial point to show that all approximations involved in deriving a theory's equations 
of motion preserve these general properties. 

One example are colloidal suspensions obeying Smoluchowski dynamics. One knows 
from general grounds that in such systems the matrices of partial correlation functions are 
completely monotone, i.e. they can be written as a superposition of decaying exponentials 
only. An approximative theory calculating such quantities should aim to reproduce these 
properties, since they are direct consequences of the structure of the time-evolution operator. 
However, the concept of complete monotonicity is quite subtle, and it is therefore likely that 
a given approximation prevents the approximative solutions from sharing this feature with 
the complete solution. 

At high densities, colloidal systems are known to undergo glassy dynamics, provided 
crystallization can be suppressed for a sufficiently long time. In these cases the so-called 
mode-coupling theory of the glass transition (MCT) has been successful in describing much of 
the experimental facts. For one-component, i.e. monodisperse, systems the theory has been 
proven [|l] to give results for the density correlation functions, which indeed reproduce the 
above mentioned features. However, little is known about multi-component, or polydisperse, 
mixtures. Recently, glassy dynamics in a binary colloidal suspension has been studied [|J, 
challenging detailed comparisons of the MCT for mixtures with experiment. 

MCT tries to describe the motion of particles on a microscopic scale by an approxima- 
tion of the frequency- and wavenumber- dependent viscosity in terms of density-fluctuation 
products. Such an approach comes about naturally if one considers the potential stresses 
to be built up by the density fluctuations of the liquid itself. The resulting equations con- 
tain a feedback mechanism induced by the slowing down of relaxations due to the so-called 
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cage effect. It has been shown that these equations allow for an ideal glass state, given the 
interactions are strong enough, e.g. at high densities. The ideal glass is characterized by a 
nonvanishing long-time limit of density correlation functions which corresponds to an elastic 
scattering contribution in the dynamical structure factor. In the "phase diagram", there 
occur critical manifolds, referred to as glass-transition singularities, that separate ergodic 
liquid from ideal glass states and can also extend into the glass state. Upon approaching 
a singularity point, the long-time limit of the density correlators, sometimes referred to as 
the glass form factor or the Edwards- Anderson parameter of the system, changes discontin- 
uously. This holds for both liquid-glass transitions, where the change is from zero to some 
nonzero value, as well as glass-glass transitions. 

Close to the transition, analytic formulas providing asymptotic solutions of the MCT 
equations of motion have been derived. There occur two time fractals with nontrivial expo- 
nents, and two diverging time scales, accompanied by corresponding scaling laws. On these 
time scales the dynamics of the glass former is also referred to as structural relaxation. The 
asymptotic predictions as well as numerical solutions of the full MCT equations for model 
systems have been extensively tested against experimental data; for a detailed discussion of 
the glass-transition scenario, the reader is referred to a recent review 0. 

The aim of this paper is to generalize from the case of one-component systems to that of 
multi-component mixtures proofs of the basic properties of MCT solutions. In particular, 
we show that the structural relaxation can be represented as a continuous superposition 
of decaying exponentials, i.e. that the density correlation functions indeed are completely 
monotone. Furthermore, the long-time limits can be obtained by a simple iteration procedure 
that does not involve the solution of the complete dynamical equations. An investigation of 
this iteration brings out glass transitions in general to be bifurcations of the cuspoid type. 
In addition, a short-time expansion of the density correlation functions is demonstrated to 
be convergent for short times, and similarly, a power series for small frequencies is shown to 
exist. 

The paper is organized as follows: In Sec. II, we will introduce the equations of motion for 
the mixtures considered, together with some basic properties and the MCT approximation. 
Section [TIT] presents a proof that the MCT equations of motion have a uniquely determined 
solution which is completely monotone. In Sec. [TV], the long-time limit of these solutions 
will be discussed, and in Sec. [V], the existence of power series solutions both in the time and 
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frequency domain is shown, the latter by proving first that all moments of the correlation 
functions are finite. Section [VT] offers some conclusions. 



II. EQUATIONS OF MOTION 
A. General Properties 

We consider a classical system enlcosed in a box of volume V with a total number of 
particles iV = X^=i ^> consisting of s different species with number concentrations x v = 
N u /N. The particles are supposed to be structureless, i.e. they do not possess any internal 
degrees of freedom, and are thus fully described by their positions, momenta, and species 
index. The variables 



are then the fluctuating densities of species v to wave vector q. Here, r% denotes the position 
of the i-th particle of species v. The simplest statistical information on structural dynamics 
that can be extracted from a multi-component supercooled liquid is the matrix of density 
correlation functions 



Here, the brackets (-|-) denote the Kubo product with (A\B) = (5 A* SB), where 5 A = 
A — (A), and (...) indicates canonical averaging. Since ^^(g,^) is the spatial Fourier 
transform of a function that is real, translational-invariant and isotropic, it is itself real and 
depends only on the magnitude of the wave vector q = \q\. The time evolution is given 
by n^(q,t) = exp[iSt]n^(q). For a liquid obeying Newtonian dynamics, the operator S is 
just the Liouville operator, which is Hermitian with respect to the Kubo product. In this 
case, time inversion symmetry implies the density correlation matrix to be symmetric with 
respect to interchange of the species indices. 

Let us focus on colloidal liquids, and S denote the Smoluchowski operator. There, time 
inversion symmetry is broken explicitly, but still the symmetry ^^(g,^) = Q utl (q,t) holds, 
although this has to be proven seperately ||. 

One then has the spectral decomposition exp[iiSt] = f exp(— 7t)<iP 7 , with eigenvalues 7 
fulfilling 7 > 0, and P 7 denoting the projector onto the corresponding linear subspace. This 




(1) 



i=i 



= (n^q,t)\n„(q)) . 



(2) 
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immediately leads to the following representation: 

$p,(q,t) = f e-^da^qi-f) , (3) 

where the measure a is concentrated on the nonnegative real axis, is symmetric in /x, u, 
and positve: da^^j) >z 0, i.e. for any set of complex numbers y u , v = l,...,n, the 
measure y* da^^j) y v is positive (summation over repeated indices is implied here and 
in the following). A function having these properties is called completely monotone. In 
particular, t) is a positive definite matrix for all times, and for all I the time derivatives 

(— l) l d\^^ v {q, t) are positive definite. The equivalence of these formulations is the result of 
the Bernstein theorem || [|. 

Let us also introduce the Laplace transform ^^(q,z) = i J °° e lz< $ M ^(g, t) dt. Then the 
representation Eq. (|3|) shows that 

®»u{q,z) = [ 1 da^j) (4) 
j z -+- 27 

is (i) analytic for z 6 C\zR~, (ii) obeys z)* = -z*), (iii) lim z ^ ioo ^(g, z) = 0, 

and (iv) Re$ M ^(g, 2) ^ for Re 2; < 0. In reverse, these four properties are enough in order 
to guarantee a representation in the form of Eq. (^) [|5], Section 5, Theorem 2.6]. The 
spectrum to) = Im z = to + iO) then is a superposition of Lorentzians 

^u(l,^) = 2 T 2 damM> ( 5 ) 
J w + 7 

which is positive as is already implied by the passivity of the system. One also has that the 
long-time limit of the correlators exists. If this quantity is nonvanishing, 

i^(g) = V(?.*- >00 )^°. ( 6 ) 

it is called the glass form factor or nonergodicity parameter, and the spectrum exhibits an 
elastic contribution 7cF flu (q)5(uj). Passivity requires F^ v (q) >z which is consistent with 
Eq. (D. 

Let us stress again that the above properties of density autocorrelation functions are 
direct consequences of the eigenvalue spectrum of the Smoluchowski operator; in general, 
they hold for any system whose time-evolution operator has nonnegative, real eigenvalues 
only. 
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B. Mode-Coupling Theory 



Mode-coupling theory starts from the formally exact representation of the density corre- 
lation matrix in terms of a memory kernel matrix. In the Laplace domain this results in the 
matrix equation 

$(g,z) = - [zS-\q) - S-\q) [ir(q) + M(q, z))" 1 ^(g)]" 1 . (7) 

Since M(q, z) — > for z — > oo, one can identify the matrices S* and r with the short time 
expansion of the time density correlation function, $(g, t) = S(q) — r~ l (q)t + 0(t 2 ). The 
matrix S(q) is called the structure factor, and from the definition Eq. (0) one checks that for 
every q it is symmetric, real and positive definite. The same properties hold for the matrix 
r(q) characterizing the initial decay of $(g). We shall throughout this paper discuss the 
above equation for q ^ and therefore assume both S(q) and r(g) to be invertible, which is 
the generic case as long as all number concentrations are nonvanishing. 

The Zwanzig-Mori formalism gives an explcit expression for the memory matrix in terms 
of so-called fluctuating forces. We shall only be concerned with its general structure, which 
in the time domain is 

M IMU (q,t) = (X^q)\1l'(t)X l/ (q))- (8) 

Explicit expressions for the fluctuating forces X^q) and the reduced resolvent 7Z'(t) can be 
worked out [[/], §. The mode-coupling approximation consists of projecting the fluctuating 
forces onto pair modes n a (k)n T (p) and factorizing four-particle correlations || [10 . 



(n a {k)n T {ft\n\t)n a ,(k')n T ,{p)) w $ CTCT ,(A?,t)$ TT ,(p,t)(5^,^ . (9) 

In order to avoid overcounting we restrict the pair modes to k > p, k' > jf with some order 
relation. In particular, for t = this implies an approximate normalization and suggests to 
introduce an approximate projector 

V 2 « Yl \n*(k)n T (p)) S^(k)S; T \(p) (n*'(k)n T ,(p)\ . (10) 

k>p 



With this projector, the MCT approximation is [|TT|, [L2| 



M™ CT (q,t) = ^{XM\nAk)n Tl {p))S-^)S^ T {p) x 

k>p 

^(^^^^^^(^^^(pXn^^Cp)!^^). (11) 
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It is useful to write this in a more transparent form by introducing super-indices s = (a, r), 

s = (o-,f): 

M^(q, t) = Y, V »M kp) Mk, t) ® $(p, t)UV V s(q, kp)* , (12) 

k>p 

where ® denotes the tensor product in the space of species indices, and the 'vertex' reads 

V^rM kp) = (Xr($\ns{k)n T ,(p))S£(k)S£(j>) . (13) 

Since the tensor product of positive definite matrices again represents a positive definite 
matrix in the corresponding product space, it is clear from Eq. (|I2] ) that M Ml/ (g, t) is positive 
definite, provided that the density correlation matrix ^^ u (q, t) is positive definite for all wave 
vectors. More generally, if we write 

F*A F * G]= l -Y^V, S [F®G + G® F) sM V*- s , (14) 
k>p 

we see that T[F, G] is symmetric in F and G, and is positive definite for every q, provided 
both F >z and G >z for every k and p. In particular, we have M(q, t) = ^ r g [ < I ) , 

Let us mention that the vertex can be evaluated and expressed in terms of the static struc- 
ture factor matrix and the three-particle static correlation functions. Usually the structure 
factor is known only approximately. Knowledge of triple correlations is often lacking entirely, 



although in principle they can be determined from computer simulation |L3]]. However, the 
property of M^ u (q,t) being a positive definite quantity is a direct consequence of the MCT 
approximation structure, and is indeed independent of the approximations made in order 
to evaluate the vertex in Eq. (0). Note in addition that one can in principle include a reg- 
ular contribution to the memory kernel, M ieg (q, z), accounting for transient dynamics not 
captured in the MCT approximation, e.g. hydrodynamic interactions and the like. Under 
the assumption of such a term being completely monotone, the following discussion remains 
valid. That there indeed exists a completely monotone solution to Eqs. (0) and (|T2"| ) shall 
be proven in the following section. 
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III. COMPLETELY MONOTONE SOLUTIONS 



A. Complete Monotonicity 

We denote the space of s x s matrices, where s is the number of species, by A. It is 
clear that A, equipped with standard matrix multiplication and Hermitean scalar product, 
is indeed a C* algebra. For elements a q G A, we form vectors a = (a q ) g= i t _M G A M . 
One easily checks that A M with all matrix operators over A defined elementwise in q, and 
equipped with the maximum norm ||a|| = max, can again be turned into a C* algebra. 
An element a G A M shall be called positive, a y 0, if a q y for every q; similarly, we use 
a y 0, or a y b, the latter meaning a — b y 0. Note that the norm preserves ordering, i.e. 
for a y b we also have ||a|| > ||6||. 

In the following, we shall assume wave vectors to be discretized to some finite set q = 
1, . . . M, such that all matrices appearing in the equations of motion are elements of A M . 

Now assume & n \t) to be completely monotone. It is then clear from Eq. (|l2[) that 
(t) inherits this property. In particular, its Laplace transform has the properties (i) to 
(iv) of Eq. (§. But then $( n+1 )(z) as defined by 

(15) 

again fulfills properties (i) to (iv). This is easily checked for (i) to (iii). Property (iv) can 
be shown in two steps. First, define K(z) by 

[ir + M in \z)] K{z) = -1. (16) 

One then has 

(ReM(z))(ReK(z)) - (t + lmM(z))(lmK(z)) = -1, (17a) 

(ReM(z))(lmK(z)) + (T + lmM(z))(ReK(z)) = 0. (17b) 

Using the second equation, one can eliminate ImK(z) in the first and find for Re z < 0, 
where ReM(z) y 0, that ReK(z) ^ 0. But we have 

[zS' 1 + S- 1 K(z)S- 1 } = -1, (18) 

and along the same lines one derives with ReK(z) ^ for Re 2; < the desired result, 
Re $(2) y 0. 
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zS 



-1 



S- 1 [iT + M {n \z)\ 1 S 



-1 



Thus, Eq. (|15D together with some completely monotone starting point, = 
exp[— (Sr)~ 1 t]S, say, defines a sequence of completely monotone functions $( n )(t), normal- 
ized to $(£ = 0) = S. Let us complete the proof of existence of a uniquely determined, 
completely monotone solution to Eqs. (|7|) and ([TJ) by showing that the thus constructed 
sequence converges uniformly to some $(£). 

B. Existence and Uniqueness 

The equation of motion in the time domain can be obtained from Eq. (0) which yields 

r$(t) + S'^it) + (M *<&)(£) = . (19) 

Here, (/ * g)(t) = f* fit — t')g{t')dt' denotes the time-domain convolution. The density 
correlation matrix is subjected to the initial condition $(£ = 0) = S, and the memory kernel 
M(t) is given in terms of the density correlators, Eq. flT2|). 

This equation of motion can be rewritten as an integral equation similar to the Picard 
equation, 

r$(t) =tS+ [ [M(t')S - S' 1 ^) - M(t')$(t - t')] dt' , (20) 
Jo 

such that the standard proof of local existence and uniqueness can be applied. In particular, 
the Picard iteration corresponding to the Laplace-domain iteration defined in the previous 
subsection is 

$("+!)(£) = S+ [ K [$ (n) (t'), $ (n) (t-f) ,& n+1 \t')] dt', (21) 
Jo 

where 

!C[x,y,z}=T- 1 (f[x}S-(S- 1 + f[y])z) . (22) 

The convergence of this iteration is proven as in the one-component case Ujfl, using a Lipschitz 
constant L for /C. If we restrict t to some finite time interval, < t < T, and the vertices V 
to some finite closed domain, we get, since ||$^(t)|| < ||5|| ensures a finite closed domain 
for x, y, z, 

\\JC(x 1 ,y 1 ,z 1 ) - JC(x 2 ,y2,z 2 )\\ < L(\\xi - x 2 \\ + \\yi -y 2 \\ + ||zi - z 2 \\) ■ (23) 
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This allows to construct a sequence X n (t) = \\^ n+1 \t) — obeying the inequal- 

ities X n (t) < 1, 

X n (t) < L /*(2X n _ 1 (t / ) + X n (t'))dt' , (24) 
Jo 

and from there, the proof of Ref. |J applies. Thus there exists for each fixed finite T a unique 
solution in the interval < t < T < oo. Furthermore, this solution depends smoothly on the 
vertices V as control parameters for any fixed finite time interval. Note that this theorem 
cannot be extended to infinite time intervals. 

Together with the preceding paragraph, we have proven the existence of a unique solution 
to Eq. (0) that it is the limit of a convergent sequence of completely monotone functions. It 
is therefore itself completely monotone, due to the continuity theorem for Laplace transforms 

HI. 



IV. GLASS FORM FACTORS 

In this section we shall prove that the glass form factor F = $(t — > oo) can be obtained 
without solving the integro-differential equation (|19|). 

Since the Laplace-transform exhibits a pole at zero frequency, &{z) = —F/z + {smooth}, 
Eq. (|7|) implies that the form factor is a solution of 

F = S-[S- 1 + N]~ 1 . (25) 

Here, N liv {q) = M^ u (q,t — > oo) denotes the long time limit of the memory kernel matrix 
which is, according to Eqs. ([i~2]) and ([i~4]) , a quadratic functional of the form factor, iV = 
T[F,F\. 

In general the coupled equations ( |25"1) and ([14]) have several solutions, e.g. F = trivially 
satisfies the equations. Since we have shown in the last section that the solution of the 
mode-coupling equation is symmetric and completely monotone, one can restrict the search 
on the positive symmetric solutions of the above equations. In this space, the solution 
corresponding to the glass form factor shall be shown to be maximal with respect to the 
semi-ordering y defined on A M . 
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A. Maximum Fixed Point 



The mode-coupling functional preserves the semi-ordering. Since F[F, G] y is satisfied 
for F, G t 0, and T is symmetric in F and G, we have N[F] —N[G] — F[F + G, F - G] ^ 0, 
if F — G >z 0. Thus one finds N[F] y N[G] if F y G. It is easy to see that inversion reverses 
the semi-ordering, i.e. for F >z G y one has G^ 1 — F -1 y 0. 

Equations ( p5|) and (|H|) suggest to introduce a continuous mapping for a set of positve 
symmetric matrices F by 

X[F] = S - [S- 1 + N[F]] ^ . (26) 

It is clear from the preceeding paragraphs that T[F] is again positive and symmetric and 
preserves the semi-ordering, T[F] y T[G\ if F y G. By induction one shows that the 
sequence F^ n+l ' = T[F^], n = 0,1,..., starting with F*- ** = S >- is monotone and 
bounded, S y F< n > h ^ (n+1) h 0, n = 1, 2, . . . , and thus converges to some fixed point 
F* y which is a solution of Eqs. and (Q). 

Suppose now there is some positive definite, symmetric fixed point F**. If we introduce 
the mapping F h-> F by F = F** + F, this maps F = F** to F = 0, and S y to 
5 = 5' — F** >- 0. The mapping is covariant in the sense that F = Z[F] holds iff F = X[F], 
provided one sets J[F] = 5 - + N[F}}~ 1 with 

iV[F] = N[F] - A^[F**] . (27) 

It is clear that the mapping X[F] exhibits the properties of I discussed above. Thus the 
sequence 

F {n) = p(n) _ p** with p(n) ag aboye 

converges to some positive definite fixed 
point F*. By continuity of all maps involved, F* = F* — F** y 0, and thus F* y F** for 
any fixed point F**. We can summarize that F* is a maximum fixed point in the sense that 
is is larger than all other positive definite, symmetric solutions of Eqs. (|25 ) and (fL4|) with 



respect to the semi-ordering introduced above. The iteration scheme defined by I converges 
to this maximum fixed point, provided the iteration is started with the upper limit S. 

B. Uniqueness and Eigenvalue 

Denote by if) the linearization of X and thus SFS. It is clear that if> is a positive linear 
map on A M in the sense that if>[f] y for all / y 0. For the one-component case considered 
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in Ref. |T[], if) corresponds to a positive matrix A qp in the sense that A qp > for all q and 
p. From the physical picture of the MCT approximation it is reasonable to assume that A 
has no invariant subspaces, and thus is an irreducible matrix. Then the Perron-Frobenius 
theorem |D| can be invoked to prove the existence of a non-degenerate, positive eigenvector 
z corresponding to the spectral radius of A. 

In the present case, the generalization of this result is guaranteed since the equivalent of 
the Perron-Frobenius theorem holds for positive linear maps on C* algebras |L6| . The phys- 
ical interpretation again leads to the assumption that if) is irreducible. Then in particular, 
the mapping if> has a non-degenerate maximum eigenvalue r, to which there corresponds a 
uniquely determined eigenvector z y 0. For all other eigenvalues A, there holds |A| < r. 
For completeness, a proof of the general Perron-Frobenius results as far as needed here is 
sketched in Appendix |A]. 

If we suppose r = 1 + S with some 5 > 0, we have, after the transformation Iwlas 
defined above, 

^z]h^z} = (l + 5)^z (28) 

with some real £ > 0. If we set = £z and define a sequence by X( n+1 ) = XfF*™)] = 
S-^FW]^- 1 forn = 1,2, . . . , we have S-SF^S z< [S' 1 + (1 + 5) F^]' 1 . But there exists 
some e > such that S' 1 + (1 + 8)F^ - S~ 1 (S~ 1 — F ( -°^)~ 1 S~ l = S' 1 + (1+8)F<® - S- 1 (S + 
SF<®S + ... )S- 1 = 5 ■ + 0(£ 2 ) y for all < £ < e, and one gets F« y F<°> y 0. 
Since 2 inherits the properties of J, it follows that S' 1 y X[F^ n+1 ^} y T[F^\ y 0. Thus 
the sequence F^ n \ n — 0, 1, . . . , is monotone and bounded, and by continuity of X converges 
to some fixed point F*. If we now choose I such that F** = F* is the maximum fixed 
point of X, y implies the existence of some fixed point y F* of X. Thus, by 
contradiction, 5 > cannot be possible, and we conclude r < 1, i.e. the maximum eigenvalue 
of if) is bounded by unity. The value of r depends on the control parameters V, and thus 
one distinguishes regular points, r < 1, from the 'critical' manifold, where r = 1. Let us 
note that F* defined as the maximum fixed point of Eq. ( p5|) exhibits bifurcation at critical 
points, identified within MCT as the ideal glass transition singularities. The non-degeneracy 
of r implies that MCT describes glass transitions in multi-component colloidal systems as 
bifurcations of the Ai type, according to the classification of Arnol'd [[17]]. This in turn 



ensures that asymptotic solutions can be worked out in the same spirit as for one-component 
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systems, with the common case being the dynamics close to a fold (A2) bifurcation 
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C. Long-Time Limit 

Now, define a dynamical mapping $ 1— ► $ similar to above by = F* + $(£). Here, F* 
shall be the maximum fixed point of 1. It is easy to see that <3? fulfills the same equations 
of motion as provided one maps T 1— > JF as discussed above. 

Since $ is a completely monotone function, the limit lim^oo $(£) = (? exists. The same 
also holds for $, and the mapping implies G = F* + G, thus G >z F*. On the other 
hand, all time-derivatives of completely monotone functions must vanish for long times, 
<9™$(t — > 00) — > 0. Therefore, one can integrate the time-domain equations of motion, 
Eq. (0), to get S^G + N[G](G - S) = 0, which is equivalent to Eq. Q2|). By this, G is a 
fixed point of Eq. (p5|), and we have G ^< F*, from which one concludes G = F* . 

Thus the maximum fixed point of Eq. (EST) corresponds to the glass form factor matrix of 
the mixture. In particular, we have explicitly generalized the iteration scheme of Ref. [|IJ that 
allows to calculate the form factors numerically without solving the full equations of motion, 
Eq. (0) or Eq. (|T9|) . Let us also note that, with the above transformation, $(t — > 00) — > 0. 

V. POWER SERIES SOLUTIONS 
A. Power Series for Short Times 

In the time-domain, there exists a power series for t < i* with some nonzero £* > 0, 



00 




(29) 



n=0 



and analogous for M(t), where m n = Ylk=o •F'l'ftk, ^n-k] ■ 
Let us suppose such a series exists. Then Eq. (|BJ) implies 



1 



n 



k\(n-k)\ 



'n+l — 



n + 1 



T 



1 



(30) 



fc=i 



for n > 1, with 0o = 5\ We prove by induction that 



O<C||0„|| < ||0o 



(31) 
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for all n. Assuming the induction hypothesis, we have for n + 1 

CI! (An I 



>n+l\ 



< 



n 



i c-ii 



"HE 



: Wn-fc P* 



fc=l 



n! 



(32) 



Since the functional JF is continuous, we can estimate ||m n || < KJ2k=o Il0fc|| H0n-fc|| < ( n + 
l)iT||0o|| with some constant Furthermore, Ylk=i k\{n — k + 1)\ / {n + 1)\ < Ylk=i V( n + 
1) < 1, and thus the proposition is proven for £* sufficiently small. 



B. Existence of Moments 



In the following, we restrict the discussion to regular points, such that the spectral radius 
of the mapping ip is less than unity, r = 1 — 2e with some e > 0, say If one writes 
SM(t)S — (ip + 5ipl§(t)})[&(t)], one has for long times, t > t with some fixed to, that Sip 
becomes arbitrarily small as $ — > 0, 5ip ^ eip. Thus one infers an upper bound for the 
memory kernel, 

SM(t)Sl(l-e)$(t)+ti(t), t>0, (33) 

with some arbitrary jj,(t) y whose Laplace transform shall be analytic for small \z\. For 
Imz > and Rez = 0, the equivalent inequality holds for the Laplace-transformed quanti- 
ties, and inserting into Eq. flj]) yields 

r< -ie$(z) r< SrS - r< SrS - ifi(z = iO) . (34) 

This shows in particular that the zeroth moment of $(£), $o — fn° ®{t)dt exists. Due to 
Eq. (H), we can, for any fixed T and 5 > 0, write j Q T ' $(t)dt r< -i(l+<J)$(-z = i/Tln(l+5)) ^ 
(1 + S)(StS — ifj,(z = i0))/e. Explicitly, we have then 

$ = StS + SM S < oo . (35) 

The existence of all other moments of $(£) and, from that, of M(t) shall be shown by 
induction. We denote these moments by $ n , 

POO 

$ n ;= / t n $(t)dt, (36) 
Jo 

and equivalently for M(t). To simplify notation, let us drop tildes in the following. 
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First, we note that any <£>(£) >z monotone and continuous such that J °° t n $(t)dt is finite 
can be written for large times as 

m = ^T, *-oo, (37) 

with some continuous K(t) >z satisfying K(t) — > for t — > oo. To see this, write for t large, 
with C some constant, ^ $(t)t n+1 = C$(t) ft /2 T n dr ^ C ^ 2 T n ^{r)dT , where the right- 
hand side vanishes as t — > oo. Since with K(t) >z 0, also ||i^(t)|| and Z(t) = max T >t ||X(T)|| 
are continuous, one can further estimate 

0< \m)\\ <|rr, *-oo. (38a) 

For n — 0, the bound can be improved. Write t$(t) r< $(t)/(lnt - lnT ) in To dr ^ 
l/(lnt — In T ) Jt^t < ^( r )^ r 5 with some T satisfying £ > T lnT > T . From this, one gets 
$(t)(f In t) ^ fT(t) In T + l^'inVo $ ( r ) rfr =< C for lon S times - 0ne concludes AT(i) ^ C/lnt, 
and specializing to T — hi t, one finds that the bound C can be made arbitrarily small, thus 
giving with obvious renaming <£>(£) = K(t)/(t\nt) with K(t — > oo) — > 0, or 

0<||$(t)||<^, f->oo, (38b) 

with continuous Z(t) ^ vanishing as t — > oo. 

Using these expressions, one infers that the convolution (M * $)(£) decays to zero faster 
than l/t n+1 for long times, 

(M *<&)(*) = |^, t^oo, (39) 

where — > for £ — > oo, given that both M and $ have the above properties and finite 
n-th moments. Write 

/•*/2 /-t/2 

||£ n+1 (M *$)(£) || < \\t n+1 / M(t')$(t - t')dt'\\ + \\t n+1 / M(£-£')$(£')dt'|| , (40) 

Jo Jo 

and both terms on the right-hand side vanish for long times. To see this, let us focus on the 

first term, 



B(t) :=t n+1 \\ ' M(t')<Z>(t-t')dt'\\ <t n+1 QnaxJ|$(T)||^ (c + / \\M(t')\\dt' 



(41) 
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with some constant C. Now for n > 1, we immediately have that \\M(t')\\dt' is bounded 
to above by some constant, and Eq. (|38a|) gives B(t) — > C ■ Z(t) —>■ for £ — > oo. Similarly, 
for n = one uses the improved bound of Eq. (|38bp together with Jj ||M(t')||<&' < C/t to 
get -B(t) — > CZ(t)/\nt (C + ln£) — > as t — > oo, which completes the proof. 

Now, suppose all moments of and M(£) up to n exist. It is clear, that then also 
the moments of (M * $)(t) up to n exist. From Eqs. ( |19D and (|3~3"D, we get £$(t) ^ u(t) — 
^[5*(M * — Sr$(t)]i which can be integrated to give, together with Eq. (j3~9| ), 

d /in+i + (n + l)S(M * $)„ - (n + l)SY$ n . (42) 

This implies the existence of both the (n + l)-th moment of and of M(t). 



C. Power Series for Small Frequencies 

Having proven the existence of all moments, we can proceed to establish a power series 
for and M(z) for small \z\. We show that the series 

i$ n (iz) n 



n=0 



111 



(43) 



converges for \z\ < zq, with zq > a non-zero radius of convergence. 

For the n-th term in the power series, Eq. fl43|) , there holds the inequality 

l $ o|| 



zV: $ I 
o II n I 



< 



n > 1 



(44) 



n! (n + l) 2 

which can be proven by induction. Taking norms on both sides of Eq. ( [42]) one obtains, 
together with || SMkS\\ < \\&k\\ + ll^fcll? 



£\\<S> 



n+l 



< 



l/^ra+ll 



[n 



+ (ll^ll + II^ID ll 5 " 1 !! H $ »-*H + ( n + ^H^ll Il $ -H • ( 45 ) 



fc=0 



But fi{z) is analytic for small \z\, so let us estimate u n < Ct^n! with some constant C > 
and some to > 0. We also choose zq such that z to < 1/2, in order to get from the induction 
hypothesis 



y n + 1 \\(b I 

-0 iTgn+ll 

(n + l)! 



< C{z t 



Zq II $q || 

(n + iy 



n / ii-j- .I 



< 5C(z t ) | 9z ||$o| 



(n + 2) 2 (n + 2) ; 



|,s- J || (2C-||<i>„||y ) - 



< 



|SV| 



gjlgolj 
{n + 2)' 



(46a) 
(46b) 
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where we have used 2~ n < 5/(n + 2) 2 and l/(n + l) 2 < 9/(n + 2) 2 and the last inequality 
holds for small enough zq. From this, the convergence of the power series of Eq. ( f43|) follows. 
Consider Cauchy sequences for m > n, \z\ < z , 



II fc | 11-2^(^ + 1)2' 

k=n k=n 

where the right-hand side becomes arbitrarily small if n becomes large. 



VI. CONCLUSION 



We have shown that, for a multi-component colloidal mixture driven by Smoluchowski dy- 
namics, the mode-coupling theory of the glass transition (MCT) provides an approximation 
to the density correlators that preserves the complete monotonicity imposed by the general 
structure of the time-evolution operator. Thus, positivity of the spectra is guaranteed. Since 
the correlation functions are superpositions of purely decaying exponentials in the sense of 
Eq. (0), the term "structural relaxation" given to the dynamics close to a glass transition 
is justified. The glass form factor was found to be determined by a maximum principle. It 
can be evaluated by an iteration scheme whose linearization fulfills the prepositions needed 
for a generalized Perron-Frobenius theorem under the natural assumption that the system 
has no decoupling wave-vector subspaces. This in turn ensures a non-degenerate maximum 
eigenvalue that was shown to be smaller than, or equal to unity for all physical states. Thus 
within the MCT of multi-component systems, the only possible glass-transition singularities 
are bifurcations of type A#, occuring at points where the Perron-Frobenius eigenvalue equals 
unity. 

The proofs presented here suggest numerical schemes for the solution of the nonlinear 
equations of motion, that are stable in the sense that they will converge to the uniquely 
defined, symmetric, completely monotone solution. In particular, the glass form factor can 
be found without solving the dynamical equations. 

Typical solutions to the MCT equations appear in Fig. [I]. We chose as the simplest case 
a binary mixture of hard spheres. Then the state of the system is determined by three 
numbers, which we take to be the total packing fraction, (p = + (ps, the diameter ratio 
S = d B /dA, and the packing contribution of the second species, v^b/v 9 - Here, the species 
are labeled with subscripts A and B, and ip a = (ndl t /6)(N a /V) are the partial volumes 
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occupied by each species. To calculate the MCT vertex, Eq. (|13D , one needs to know the 
static structure factor of the system, as well as static three-particle correlation functions. 
For the latter, to our knowledge, no analytic expressions are available, and thus we follow the 
commonly applied approximation through two-particle correlations [IT]. The former shall 
be approximated using the well-known Percus-Yevick (PY) approximation to the Ornstein- 
Zernike integral equation. Within this framework, "PY-exact" solutions are available |I9 



As was already mentioned in connection with Eq. (|l3l), the approximations involved in 
evaluating the vertex are of no importance for the mathematical aspects of the solutions we 
wish to demonstrate. 

Figure [l] shows normalized correlation functions (j)(q,t) = S~ 1 / 2 (q)&(q,t)S~ 1 / 2 (q) for a 
particular wave vector q = 6.2/ d&. They are numerical solutions of Eq. (|7]) with T a p{q) = 
l/{q 2 D° a )5 a p and D° a = 0.005/ R a , on a grid of 100 wave vectors q = 0.2, . . . 39.8 with Aq = 
0.4. The size ratio and the concentrations are kept fixed at (8, fs/'-P) — (0.6,0.2), and the 
total packing fraction is varied as indicated in the figure caption; a procedure corresponding 
to what has been done in experiment [[J . Note that while the two diagonal matrix elements 
are positive and monotonically decreasing for all times t, both needs not be true in general 
for the off-diagonal elements. For comparison also shown is the suggested starting point of 
Eq. (|T5|), = S~ l l 2 exp[— (Sr)~ 1 t]S 1 / 2 , which is completely monotone by construction. 

One notices a drastic slowing down in the relaxation of the correlation functions towards 
their long-time limits. This demonstrates the dynamics close to a critical point where the 
Perron- Frobenius eigenvalue r = 1, which in this system occurs at tp c m 0.519608 and which 
is the concern of the asymptotic solutions to mode-coupling theory. The solutions for the 
long-time limit at tp = tp c as determined from Eq. (]25|) are included in the figure. 

Complete monotonicity requires all eigenvalues of <p(t) to be positive for all t. To demon- 
strate that our numerical solution is in accordance with this, we show in Fig. |2] the eigenvalues 
e±(t) as functions of time for a fixed total packing fraction. One clearly recognizes the above 
statement to hold. 

A power series in the time domain exists, and at regular points, i.e. for vertices such that 
the Perron- Frobenius eigenvalue is smaller than unity, also for small frequencies the power 
series has a nonzero radius of convergence. In particular, the existence of all moments of 
the density correlation function at regular points, and of a finite radius of convergence for 
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the frequency-domain power series implies the existence of a final exponential relaxation, 

$(t) — F* — C(e- 70 ') . (48) 



This holds, since Eq. ( f4"3|) implies that the measure a in Eq. ([3]) has an atom of mass F* at 
7 = 0, is constant for < 7 < 70, and has a point of increase at 7 = 70 [13]. Thus, 

/■oo 

«W9. = F*,M + / ^ do Wfl ( 7 ) • (49) 

■'70 

The mode-coupling theory for mixtures can also be applied to systems with Newtonian, 
instead of stochastic, short-time dynamics, e.g. metallic melts. Furthermore, a recent ex- 
tension of MCT to molecular liquids that treats each molecule as consisting of s constituent 
sites also leads to substantially the same equations [p0| . Since in these cases, the represen- 
tation of Eq. (|3]) through decaying exponentials only will not be valid in general, the proofs 
presented here cannot readily be applied. Work published for one-component Newtonian 
systems [[H] suggests that existence and uniqueness of the solution, even if it will not be 



completely monotone, can nevertheless be proven. This has, however, not been done so far. 

It is an observation of both theory [^] and computer experiment |23|] , that the different 
short-time behavior does not influence the dynamics at sufficiently long times apart from an 
overall shift in time scale, given a strong enough coupling such that MCT contributions are 
important. One then expects the long-time limit of the Newtonian dynamics solutions to 
exist and to be governed by Eq. (|25|) . Note that the properties of this equation and its max- 
imum fixed point do not depend on the short-time dynamics, nor do the commonly applied 
asymptotic formulas. Similarly, the prediction of only Ag singularities as glass transitions 
will remain valid as long as the linearization of Eq. (|25|) is irreducible. This can be expected 
unless some special symmetry will introduce zero-couplings in the vertex, which in principle 
can happen within the molecular site-site description of Ref. [2(J . 
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APPENDIX A: PERRON THEOREM 



Let us sketch here for completeness some results generalized from the Perron- Frobenius 



theorem for irreducible matrices ||15|| . For a generalization of the complete theorem for 
positive linear maps on C* algebras, the reader is referred to Ref. ||16|| . 

As above, let A M denote the C* algebra of M-vectors of s x s matrices over C. Consider 
the positive linear map ip, which maps the set of symmetric, real, positive definite elements 
At. onto itself, ip[ a } ^ for a >z 0. ip is called irreducible, if there exists some positive, finite 
number n such that T[a] y for a y and 

T[a] = (l + ^)>]. (Al) 

If i/j is irreducible, we have that ip[a] y if a y 0. 
Now, define a mapping r : (A¥ , C s ) — > K by 

r a,v = min v ' , gL \ ' , A2 
i< 9 <M (v\a q v) 

where q labels the elements of a e ^4 M , a = (a q ) q =i,...M, and (-|-) is a scalar product over C s . 
Furthermore, set r(a) = inf^gc^ r(a,t>) = inf„ e§s r(a, v), where S s denotes the s-dimensional 
unit sphere, and the latter equation holds since r(a, Xv) with A G C is independent of A. 
One immediately checks r(a) > and 

^]br(o)o. (A3) 

However, r(a,v) is not continuous on (A^,C S ). Let us define a set 

B := {6; b = T[a},a G ^Jf , ||a|| = 1} . (A4) 

Then, 6^0 for any b G B C ^4^ 7 . Since r(b,v) is continuous on the closed and compact 
set (£>, it assumes its minimum with respect to v. It follows that on £>, r(6) fulfills a 
maximum principle in the sense that it is the maximum IR-number for which ip[b] >z r{b)b. 
Furthermore, for a G A+ and b = T[a] G B we have T[ijj(a) — r(a)a] = i/}[b] — r(a)b y 0, and 
by the maximum principle we get r(b) > r(a). 
Next define 

r = sup r(a) . (A5) 
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Since r(l) > 0, clearly r > 0. The supremum can be restricted to b E B since the above 
inequality holds. But there, r(b) assumes its maximum, and thus r attains the supremum 
for some extremal vector z >- 0. 

We continue by showing that r indeed is an eigenvalue of ip and equal to the spectral ra- 
dius, and that the corresponding eigenvector z is unique, i.e. the eigenvalue is non-degenerate. 

Assume, ip[z] — rz >z but not the null element. Then ip[z] — rz >- with z, = T[z], 
but the maximum principle then implies r(z) > r in contradiction to the definition of r. 
Thus, r is an eigenvalue of tp. Suppose now, there are two eigenvectors z, z' corresponding 
to this eigenvalue. We then can find some A G R such that \z — z' >z but not strictly 
positive definite. But this implies T[Xz — z'\ = (1 +r) n (\z — z') >- 0, in contradiction to the 
construction of A. Thus, the eigenvalue r is non-degenerate. 

Now for any a G A+ , define the mapping 

a[a] = (l/^z-^lz^az 1 ^- 1 / 2 . (A6) 

Since a[l] = 1, \\a\\ = 1. Suppose ip[u] = au with some a G C. Write v = z~ x l 2 uz~ 1 l' 1 
which gives a[v] = (a/r)v. But for any eigenvalue A of a, we have |A| < \\cr\\, and thus in 
particular |a| < r. Therefore, r is the spectral radius of ip. 
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FIG. 1: Normalized correlation functions (t> a p(q,t) = (5 _1 / 2 (g) < I ) (g, t)S~ l l 2 {qj) a p for a binary 
hard-sphere mixture with diameter ratio d-^/dx = 0.6, packing fraction of B-particles <Pb/<P = 0.2, 
wave vector q = 6.2/dA, calculated from Eq. ([?]); details see text. The labels 11, 12, and 22 
indicate the three independent matrix elements, and labels a and b correspond to total packing 
fractions ip = 0.5195 and 0.525, respectively. Here, (p = {-kN /QV){xxd\ + xb^b), with x a = N a /N 
being the number concentrations. The dashed horizontal lines indicate the long-time limit for 
(p = (p c « 0.519608, and the dotted curves correspond to an exponential decay as described in the 
text. 
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FIG. 2: Eigenvalues e±(q,t) of (j>{q,t) from Fig. [l] for tp = 0.5195 (all other parameters as given in 
Fig. [I]) as functions of time t. 
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